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O ' Abstract 

£>y The optimization of the light scattered by photonic cluster made 

« . of small particles is studied with the help of the local perturbation 

method and special optimization algorithm. It was shown that pho- 
tonic cluster can be optimized in a such a way that its reflectivity will 
be increased or decreased by several orders of magnitude for selected 
^_^ ' wavelength and direction. 

in ' 
in 

)Q 1 Introduction 

o 

Light is known as fastest carrier of energy and information and this prop- 
erty makes light indispensable for communications, warfare and fundamental 
research. While light can be relatively easy created and guided, its manip- 
ulation is somewhat difficult. Manipulation of the light requires dynamical 
control of the refractive index of the host medium. While the materials for 
active control of the light are still under development, the theoretical stud- 
ies are already started. Recently, the light manipulation was investigated in 
work m where the weak scattering by the finite object was studied by using 
the first Born approximation. In practice, the light scattering by photonic 
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cluster may be not weak and it should be studied with other methods. The 
local perturbation method (LPM) is suitable tool for such analysis. The LPM 
correctly describes scattering by particles with arbitrary large refractive in- 
dex which are small compared to the incident wavelength and the scattering 
can be strong (see for example work [2]- [4] and references wherein). In works 
[5]- [7] the LPM was used to study the wave propagation in the photonic 
cluster. 

In this paper we optimize the light scattering from the cluster made of 
small particles by using the LPM and special optimization technic. By using 
our method we modify the cluster in such a way that the scattering from 
the cluster is significantly minimized or maximized at one point. We present 
several examples which will show that the scattering by the cluster can be in- 
creased and decreased by several orders of magnitude for selected wavelength 
and direction. 

2 The formalism 

The formalism we use is presented in many works (see for example [2]- [I]) 
and we only briefly present it here for convenience and consistency. Consider 
the photonic cluster made of particles which characteristic sizes are small 
compared to the incident wavelength Ao- The electric field E propagating in 
the host medium filled with N particles is described by the following equation 

® 



o JV-1 



c^ 



(A - V <g> V + k 2 ) E(r) + — J> sc , n - e„)/„(r - r n )E(r n ) = S(r), (1) 
where 



n=0 



2-7T uj — , . ( 1, inside particle 

A c ' \ 0, outside particle 

Here A is Laplacian and V nabla operators, <8> defines tensor product, fco 
is a wave number in the host medium (u is the angular frequency and c is 



the light velocity in vacuum), e sc>n and Eq are the permittivity of the n-th 
particle and the medium respectively, f n is the function describing the shape 
of the n-th scatterer, and S is the field source. The characteristic size of n-th 
scatterer we denote as L n . Note, that the equation flTJ is an approximate one 
and this is true only when the small scatterers (k L n <C 1) are considered. 
This equation is easily solvable with respect to the fields E(r n ) when the 
positions r n are known (it will invoke solution of 3N linear equations for 
each frequency uj). The solution of the equation ([[]) can be written in the 
following form 

E(r) = E (r) + E sc (r), (2) 

where 

E sc (r) = ^ 1 1 + Z®Z j J2 E(r n )(e sc>n - e )<S> n (r) (3) 

V o / n=0 

and 



*" (r) = L (2-*8) " q - Uq> = S? /„ / " (r) ^"" dr ' (4) 

Here J is the 3x3 unitary tensor in polarization space and r n is the radius 
vector of the n-th particle. The field E(r n ) is the field inside the n-th particle, 
/„ is the Fourier transform of the function f n , and • defines scalar product. 
The incident field E is created by the source S in the host medium and it 
is not important for our consideration. 

The formula (T5J) is rather general one and it describes the field in the 
medium with photonic cluster of arbitrary form made of small particles of 
arbitrary form. 

The fields E(r n ) should be found by solving the system of 3iV linear 
equations obtained by substituting r = r n into Eq. (|2J). The formula for the 
scattered field ([3]) can be simplified even further when the distance between 
the observer and an n-th scatterer (R n ) is large, i.e. when R n ^> L n . In 
this case the integral $ n can be calculated approximately. We note also that 



the integral $ n can be calculated exactly at least for the spherical particles. 
When R n 3> L n integration in (J3J) gives 

^{r) = -^[l+-^-) "£v(r n )(e sc , n -e )—V n , (5) 



where 

R n = \r- r n | > L n . (6) 

Here i? n is the distance between the observation point and n-th scatterer, V n 
is the volume of the n-th scatterer, and |...| brackets denote absolute value. 
In practice, the distance between the cluster and the observer is much 
more larger than the size of the cluster and the inequality |r| ^> max(|r n |) is 
fulfilled. In this case the field (jSJ) can be simplified and rewritten in the form 



E sc (r) 
where 
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(T - 1 <g> l) Y, E(r n )V n (e«, n - e )e^ k ^\ (7) 



n=0 



1 = r/r, r = |r| ^> max(|r n |). 

3 Optimization of the cluster reflectivity 

The idea of our approach is the following. We assume that scatterers in the 
cluster can be repositioned in a such a way that the field scattered by the 
cluster will be tuned to required value at the observation point. 

At first, we calculate the intensity of the scattered field from the given 
(not modified) cluster for selected wavelength \ opt and observation point r. 
After this we select one particle in the cluster and reposition it a bit in x, 
y, and z directions (usually for tenth of the particle size) recalculating the 
scattered field in r each time we move the particle. From the array of the 
calculated intensities we select the value nearest to the required one and we 
place the particle at the point corresponding to this intensity value. We 
perform this procedure for all particles in the cluster and repeat it several 
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Figure 1: The logarithm of the normalized intensity of the scattered field 
J sc (r)//j nc (0) versus k L for the not optimized Au cluster. The particles 
in the cluster are arranged into cubic lattice. The permittivity of the host 
medium is Eq = 1, the characteristic size of the particles is L = 10 nm, the 
period of the cluster is d = 3.3L, and the total number of the particles is 
N = 123. The source is positioned at r s = {a; s ,0,0} and the observation 
points r are located at the circle with radius x a . The angle 9 is the angle 
between vector r and r s in xy plane. 



times. After several iterations, the particles in the cluster will be rearranged 
in such a way that the intensity of the scattered light will correspond the 
required value at the selected observation point r and for the wavelength 

A-opt- 

For example, super reflective cluster can be created when the intensity of 
the scattered field I sc = |E sc (r)| will be maximized at some point and min- 
imized at all other points (directions). Another important example closely 
related to optical cloaking is the cluster with significantly reduced reflectivity 
for selected wavelength at chosen direction. When the reflectivity of the clus- 
ter is lower than the sensitivity level of the receiver, this cluster is actually 
invisible for observer. 




Figure 2: The logarithm of the normalized intensity of the scattered field 
I sc (r)/I inc (0) versus k L for the optimized spherical cluster made of Au cubes. 
The permittivity of the host medium is Eq = 1, the characteristic size of the 
cubes is L = 10 nm and the total number of the particles is N = 123. 
The source is positioned at r s = {x s ,0,0} and the observation points r are 
located at the circle with radius x s . The angle 9 is the angle between vector 
r and r s in xy plane. The cluster reflectivity was minimized for k§L = 0.1 
at r = r s . The red line shows extremely low deep in cluster reflectivity after 
the optimization. 



Moreover, when the intensity of the field scattered by the cluster is 
changed in some direction, this can be used in directional optical switch. 
This kind of device can be used for light houses, for example. When the 
spectrum of the field scattered by the cluster is modified in time when the 
cluster can be used as a filter. 

Below we present several examples demonstrating the possibility to opti- 
mize the cluster reflectivity in chosen direction for selected wavelength. 

We note that similar approach was successfully used in work [8] where 
the phase of incident field was tuned to maximize the field scattered from 
the complex object (eggshell and TiO powder) in chosen direction. 



3.1 Examples 

Consider the spherical photonic cluster consisting of Au cubes. The charac- 
teristic size of the cubes is L and they are organized into simple cubic lattice 
with period d. The center of the cluster is positioned at the beginning of the 
coordinates. The incident field is generated by the point source positioned 
at the point r s = {x s ,0, 0} and it is linearly polarized in y direction. The 
scattering by the cluster will be optimized (to minimum or maximum) for the 
chosen wavelength X opt for the observer positioned at the point r. We assume 
that the permittivity of the scatterers is the same as the permittivity of the 
bulk Au and the actual values of permittivity were taken from [9] . Note that 
the source and the observer are far from the cluster such that k \r s \ ^> 1 and 

fco|r| > 1. 

The examples of the optimization are presented on Figs. [T]l5j The figures 
show the normalized intensity of the light scattered by the clusters versus 
koL. The Fig. [T] shows the light intensity scattered from the not optimized 
cluster. Figs. [2] and [3] show the results for the clusters which scattering was 
optimized at the point r = r s and the Figs. H] and [5] show the results for the 
clusters which scattering was optimized at the point r = {0,x s ,0}. 

The positions of the particles in the clusters were changed in order to 
maximize or minimize the scattering for the wavelength X opt = 633 nm 
(koL = 0.1) in chosen directions. The results are presented for different 
angles of observation 6 in xy plane. The Fig. [2] shows normalized (with 
respect to the intensity of the incident field li nc {0)) intensity of the light 
I sc (r) = |E sc (r)| scattered by the optimized photonic cluster. The cluster 
was optimized to minimal intensity of the scattered field in the direction of 
the source (6 = 0). Comparing Fig. [TJ and Fig. [2] one can see huge deep 
(about five orders of magnitude) in the intensity of the scattered field for 
the selected wavelength X opt and the selected direction 9 = 0. The Fig. [2] 
shows also that unintendently, the optimization significantly increased inten- 
sity of the scattered field for other wavelengths (near k L = 0.07) in almost 
all directions. 

The Fig. [3] shows normalized intensity of the light scattered by the cluster 
optimized to maximal intensity of the scattered field in the direction 6 = 0. 
Comparing Fig. [T]and Fig. [3] one can see relatively large increase (about two 
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Figure 3: The logarithm of the normalized intensity of the scattered field 
/ sc (r)// inc (0) versus k L for the optimized spherical cluster made of Au cubes. 
The permittivity of the host medium is e$ = 1, the characteristic size of the 
cubes is L = 10 nm and the total number of the particles is N = 123. The 
source is positioned at r s = {x s , 0, 0} and the observation points r are located 
at the circle with radius x s . The angle 9 is the angle between vector r and 
r s in xy plane. The cluster reflectivity was maximized for k L = 0.1 at 
r = r s . The red line shows extremely high peak in cluster reflectivity after 
the optimization. 



times) in the intensity of the scattered field for the selected wavelength \ opt 
and the selected direction 6 = 0. The Fig. [3] shows also that the optimization 
significantly increased intensity of the scattered field for k L = 0.09 and 
significantly decreased the intensity for k L = 0.11 in directions with 9 = n/2 
and 6 = 7r/3 respectively 

The Fig. H] shows normalized intensity of the light scattered by the cluster 
optimized to minimal intensity of the scattered field in the direction perpen- 
dicular to the source-cluster direction (6 = ir/2). Comparing Fig. [T]and Fig. 
H] one can see significant decrease (about two orders of magnitude) in the 
intensity of the scattered field for the selected wavelength X opt in the selected 
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Figure 4: The logarithm of the normalized intensity of the scattered field 
l sc (r)/ij nc (0) versus k L for the optimized spherical cluster made of Au cubes. 
The permittivity of the host medium is Eq = 1, the characteristic size of the 
cubes is L = 10 nm and the total number of the particles is N = 123. 
The source is positioned at r s = {x s ,0,0} and the observation points r are 
located at the circle with radius x s . The angle 9 is the angle between vector 
r and r s in xy plane. The cluster reflectivity was minimized for k§L = 0.1 at 
r = {0,x s ,0}. The red line shows extremely low deep in cluster reflectivity 
after the optimization. 

direction 9 = ir/2. In distinction to previous examples, the intensity of the 
scattered field in other directions is almost not affected in this case. 

The Fig. [5] shows normalized intensity of the light scattered by the cluster 
optimized to maximal intensity of the scattered field in the direction 9 = 7r/2. 
Comparing Fig. [TJand Fig. |5]one can see huge increase (about three orders of 
magnitude) in the intensity of the scattered field for the selected wavelength 
X op t in the selected direction 9 = n/2. We note that the intensity of the 
scattered field increased also in direction 9 = 7r/3 (for A opt ) and for all other 
directions the intensity peak shifted to k L = 0.09. 

Examination of the Figs. [T]l5] reveals several important things. The first 
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Figure 5: The logarithm of the normalized intensity of the scattered field 
I sc (r)/I inc (0) versus k L for the optimized spherical cluster made of Au cubes. 
The permittivity of the host medium is Eq = 1, the characteristic size of the 
cubes is L = 10 nm and the total number of the particles is N = 123. 
The source is positioned at r s = {x s ,0,0} and the observation points r are 
located at the circle with radius x s . The angle 9 is the angle between vector 
r and r s in xy plane. The cluster reflectivity was maximized for k$L = 0.1 at 
r = {0, x s , 0}. The black line shows extremely high peak in cluster reflectivity 
after the optimization. 

one is that the optimization works because the intensity of the light scat- 
tered by the optimized cluster has very clear minima or maxima for the 
chosen wavelengths and positions. The second one is that the intensity of 
the scattered field has deep or peak not only for designed wavelength but for 
other wavelengths and directions. It should be emphasized that the scatter- 
ing was optimized for one direction and for one wavelength only and that 
is why the optimization is effective for narrow spectrum of wavelengths and 
scattering directions. 
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4 Conclusions 

The reflectivity of the spherical photonic cluster made of small particles was 
optimized to maximal and minimal values by using the local perturbation 
method and the special optimization algorithm. The possibility to design the 
photonic cluster with required scattering characteristics was demonstrated 
for selected wavelength and direction. 
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